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(57) Abstract: The preferred embodiment provides a system and method for geo referencing digital raster maps by associating 
specific points on the raster map with corresponding points on a georeferenced vector map or another georeferenced raster map. 
According to the preferred embodiment, a raster map and a corresponding vector map are simultaneously displayed to a user. The 
user then locates a common geographic point or feature on each map, and marks each of them as a unique point-pair. When the 
user has marked at least two point-pairs, the system then computes a georeferencing function, based on the pixel-coordinates of the 
points marked on the raster map and the corresponding geographic coordinates of the points on the vector map. Thereafter the geo- 
graphic coordinates of any point on the map may be easily computed. The preferred embodiment provides that as more point-pairs 
are identified the tjeoreferencine function is modified for inc-eased accij-^rv 
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SYSTEM AND METHOD FOR GEOREFERENCING DIGITAL 

RASTER MAPS 

CROSS REFERENCE TO RELATED APPLICATIONS: 

This application contains some specification and figures in common 
5 with concurrently filed, copending applications "System and Method for 

Performing Flood Zone Certifications'* (09/537,161) filed March 29, 2000 and 
•*S>stcm and Method for Synchronizing Raster And Vector Map Images" 
(00 537.162) filed March 29, 2000 which are hereby incorporated by 
reference. 

10 BACKGROITND OF THE INVENTION 

Tecbnieal Held: 

The present invention generally relates to geographic map processing 
and in particular lo a system and method for associating digital raster maps 
with known geographic features. Still more particularly, the present invention 
15 relates to a system and method for associating specific points on digital raster 

maps u iih a geographic coordinate system. 

Description or the Related Art 

A digital raster map is a computerized map image that resembles a 
conventional paper map in that it presents an image of the mapped area, but 

20 has no additional underlying data associated with the features of the map. A 

raster map is t>pically created by scanning a conventional paper map, and is a 
grid-based map composed of pixels (or dots) of color or black & white. Each 
pixel in the grid has can be referenced by the pixel coordinates, and has only 
one value, which indicates the color of that pixel. Raster images are 

25 commonly referred to as "bit mapped." 

A vector map uses lines and polygons, rather than pixels, to describe 
an image. Unlike a raster map, in which the map image is simply stored as a 
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grid of pixels, when a vector map is displayed, it is drawn from a set of 
underlying data. The vector map is created through the placements of nodes 
on a plane and connecting those nodes with lines. Vector lines can be 
attributed with tables of data such as elevations, values, names or other 
5 information relative to the line. Vector data can be displayed in three 

dimensions if the lines are attributed with z values, modified or changed 
relative to user need, or layered to allow for turning off and on the viewing of 
different information. 

Because of their feature attribution properties, they are particularly 
10 useful for displaying geographic data. Vector maps are used to display 

boundaries or lines that denote the position and extent of features, such as 
county boundaries or lines denoting stream and river systems. It is also very 
easy to view or manipulate the data underlying a vector map, for example to 
view or change the elevation of a feature. 

15 Also because vector maps are commonly drawn from a geographic 

data set describing the area shown, they are very easily, and generally 
inherently, georeferenced. Georeferencing is the process of relating source 
coordinates to referenced geographic coordinates., which are typically in 
standard latitude/longitude. An image or a vector file is georeferenced to be 

20 used within a mapping/geographic environment. In a vector map, the data 

from which the map is drawn will typically already include a geographic 
coordinate set. 

Modem GIS systems normally make use of digital vector based map 
information. However, a vast legacy of paper-based map information exists. It 
25 is very expensive and time consuming to convert all of the information on 

these paper maps over to a digital vector format. In many cases the scope and 
expense of such conversions render them completely impractical. However, 
even when a complete conversion to digital vector-based format is not 
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possible, it is still possible to obtain some of the benefits of computerized map 
systems, first by converting the paper maps to digital raster maps (by scarming 
them), and then by georeferencing the raster image. After georeferencing, 
there should be a clear relationship between the pixel coordinates in the raster 
5 map, and the geographic coordinates of the feature represented by that pixel. 

It would be desirable to provide a system and method for georeferencing a 
raster map by associating points on that map with corresponding points on a 
previously-georeferenced vector map, or with reference points designated in 
the raster map which have known latitude and longitude. 
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SUMMARY OF THE INVENTION . 

It is therefore one object of the present mvention to provide improved 
geographic map processing. 

It is another object of the present invention to provide an improved 
s\T>tcm and method for associating digital raster maps with known geographic 
features. 

It is yet another object of the present invention to provide an improved 
s>'stcm and method for associating specific points on digital raster maps with a 
geographic coordinate system. 

The foregoing objects are achieved as is now described. The preferred 
embodiment provides a system and method for georeferencing digital raster 
maps by associating specific points on the raster map with corresponding 
points on a georcferenced vector map or another georeferenced raster map. 
According to the preferred embodiment, a raster map and a corresponding 
vector map are simultaneously displayed to a user. The user then locates a 
common geographic point or feature on each map, and marks each of them as 
a unique point-pair. When the user has marked at least two point-pairs, the 
system then computes a georeferencing function, based on the pixel- 
coordinates of the points marked on the raster map and the corresponding 
geographic coordinates of the points on the vector map. Thereafter the 
geographic coordinates of any point on the map may be easily computed. The 
preferred embodiment provides that as more point-pairs are identified, the 
georeferencing function is modified for increased accuracy. 

The above as well as additional objectives, features, and advantages of 
the present invention will become apparent in the following detailed written 
description. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features believed characteristic of the invention are set forth 
in the appended claims. The invention itself however, as well as a preferred 
mode of use, further objects and advantages thereof, will best be understood 
5 by reference to the following detailed description of an illustrative 

embodiment when read in conjunction with the accompanying drawings, 
wherein: 

Figure 1 depicts a data processing system in accordance with a 
preferred embodiment of the present invention; 

10 Figure 2 is an image of an exemplary raster map, in accordance with 

the preferred embodiment; 

Figure 3 is an image of an exemplary vector map, corresponding to the 
raster map of Figure 2, in accordance with a preferred embodiment of the 
present invention; and 

15 Figure 4 is a flowchart of a process in accordance with a preferred 

embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

With reference now to the figures, and in particular with reference to 
Figure 1, a block diagram of a data processing system in which a preferred 
embodiment of the present invention may be implemented is depicted. Data 
5 processing system 100 includes processor 102 and associated L2 Cache 104, 

which in the exemplary embodiment is connected in turn to a system bus 106. 
System memory 108 is connected to system bus 106, and may be read firoin 
and written to by processor 102. 

Also connected to system bus 106 is I/O bus bridge 110, In the 
10 exemplary embodiment, data processing system 100 includes graphics adapter 

118 coimected to bus 106, receiving user interface information for display 
120. Peripheral devices such as nonvolatile storage 114, which may be a hard 
disk drive, and keyboard/pointing device 116, which may include a 
conventional mouse, a trackball, or the like, are connected to I/O bus 112. 

15 The exemplary embodiment shown in Figure 1 is provided solely for 

the purposes of explaining the invention and those skilled in the art will 
recognize that numerous variations are possible, both in form and function. 
For instance, data processing system 100 might also include a compact disk 
read only memory (CD-ROM) or digital video disk (DVD) drive, a sound card 

20 and audio speakers, and numerous other optional components. All such 

variations are believed to be within the spirit and scope of the present 
invention. Data processing system 100 is provided solely as an example for 
the purposes of explanation and is not intended to imply architectural 
limitations. 

25 The preferred embodiment provides a system and method for 

georeferencing digital raster maps by associating specific points on the raster 
map with corresponding points on a georeferenced vector map or another 
georeferenced raster map. According to the preferred embodiment, a raster 
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map and a corresponding vector map are simultaneously displayed to a user. 
The user then locates a corrunon geographic point or feature on each map, and 
marks each of them as a unique point pair. When the user has marked at least 
two point-pairs, the system then computes a georeferencing function, based on 
5 the pixel-coordinates of the points marked on the raster map and the 

corresponding geographic coordinates of the points on the vector map. 
Thereafter the geographic coordinates of any point on the map may be easily 
computed. The preferred embodiment provides that as more point-pairs are 
identified, the georeferencing function is modified for increased accuracy. 

10 Figure 2 is an exemplary raster map, in accordance with the preferred 

embodiment. This exemplary map shows a scanned image fi-om a Federal 
Emergency Management Agency (FEMA) paper map. This raster image 
shows a land area with flood zone indications, but would, in a computer 
system, contain no underlying data regarding the area shown. 

15 Figure 3 is an exemplary vector map, corresponding to the raster map 

of Figure 2, in accordance with a preferred embodiment of the present 
invention. This map shows the same area as the map in Figure 2, but is 
created by a computer system from a database describing the locations of 
features such as the streets shown. Typically, each feature shown on a vector 

20 map such as this will already be georeferenced, in that the geographic 

coordinates of each feature will also be recorded in the underlying data 

A digital map image is considered georeferenced if a pair of 
mathematical functions, f, and g, have been defined that can be used to convert 
back and forth between the coordinates of the map image (as defined by the 
25 pixels of the image) and the corresponding longitude and latitude of the 

location of that point. That is, f and g do the following: 

1 . If (x, y) represents a location on the digital map image, then 
f(x, y) =(Lon, Lat) represents the longitude and latitude of the corresponding 
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physical location. 

2. If (Lon, Lat) represents a physical location that lies within the 
region covered by the map, then g (Lon, Lat) = (x, y) represents the point on 
the digital map image that corresponds to that longitude and latitude. 

5 Here, x and .y represent the natural internal coordinate system of the 

map image. Typically, as described above, a digital raster map image uses the 
pixels of its image as a natural coordinate matrix. However, in most cases, a 
vector-based map image uses longitude and latitude as its internal coordinate 
system; if so, it can be considered to be trivially georeferenced already. 

10 Therefore, the functions f ( ) and g ( ) above are non-trivial georeferencing 

fimctions required to convert back and forth between coordinate systems. 

Following is a description of a georeferencing process in accordance 
with the preferred embodiment: 

First, the raster map to be georeferenced (referred to as Mapl) is 
15 displayed on a computer monitor in conjimction with another previously 

georeferenced map (referred to as Map2) (step 400). Typically, Map2 will be 
a vector map, but could also be another raster map, for which georeferencing 
functions are already known. Mapl and Mapl can be independently 
manipulated for rotation, position and scale. The user will manipulate the 
20 maps until the region shown in Mapl is contained within the region shown by 

Map2. 

Next, pairs of points representing identical geographic locations on 
Mapl and Map2 are identified, and the user uses a keyboard, mouse, or other 
pointing device to mark these point-pairs (step 405). For each point-pair, one 
25 point is marked on Mapl and the other point is marked at the corresponding 

location on Map2. Examples of common georeferencing point-pairs include 
street intersections, places where rivers cross roads, mountain peaks, building 
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locations, and other readily identifiable map features. If any points on Mapl 
have points of known longitude and latitude already marked on the map, then 
a georeferencing point-pair can be obtained by marking that point on Mapl and 
moving to the point with identical latitude and longitude on Map2 and 
5 marking it as the other half of the point-pair. Alternatively, once the known 

points are marked with the known latitude and longitude, they can be used to 
deteraiine the georeferencing functions without requiring a corresponding 
markonMap2. 

After a point-pair is marked, the x-y point on Mapl is assigned the 
10 corresponding latitude and longitude of its matching point on Map2 (step 

410). 

When two georeferencing point-pairs are marked, the system will 
compute georeference functions for Mapl based on a linear transformation 
that allows an arbitrary rotation and predefined scaling, between the x-y 
15 domain and the latitude- longitude domain (step 415). 

To improve accuracy and allow increased generality, more 
georeferencing point-pairs will firequently be desired- To facilitate the process 
of gathering these additional georeferencing point-pairs, Mapl and Map2 are 
now selectively synchronized. This means that as either map is zoomed, 

20 panned, scrolled, or otherwise caused to display a different region, then the 

other map automatically does the same. This greatly increases the ease with 
which georeferencing point-pairs can be identified and marked on the maps, 
since when the user of the system locates an identifiable feature on one map, 
the other map will automatically be showing very nearly the same location on 

25 the other map. 

Thereafter, when the user marks a georeferencing^oint on one of the 
maps, then the system automatically marks the corresponding point on the 
other map, based on its current georeferencing function (step 420). The user 
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can then make corrections to the location of a point marked by the system in 
this way (step 425), The user may enter as many point-pairs as are desired; 
each additional point-pair will increase the accuracy of the georeferencing 
process. 

5 When three or more georeferencing point-pairs are determined, the 

system computes georeference functions for Mapl based on a completely 
general linear transformation between the x-y domain and the latitude- 
longitude domain. 

When four or more georeferencing point-pairs are determined, the 
general linear georeferencing functions are over-determined. This meaiis that 
more than the required amount of information to compute the general linear 
georeferencing functions is available, but that it is not, in general, completely 
consistent. The system use the extra information contained in the additional 
georeferencing points to provide validation checks to protect against the 
possibility that some of the data points may be inaccurate (step 430). Points 
that deviate excessively with respect to a calculated standard error are 
presumed to be inaccurate and are omitted from the calculation of the 
georeferencing functions. Note that as new point are added, the system also 
rechecks points previously marked as inconsistent, to determine if those 
points should now be considered when recomputing the georeferencing 
functions. 

The additional point-pairs allow the system to compute the general 
linear georeferencing functions which best fit the combination of all the 
available georeferencing point-pairs, so that the active georeferencing 
25 functions are revised with each new point-pair (step 435). According to the 

preferred embodiment, this is done by using a "least square" parameter fitting 
operation. 

The user may then proceed o enter the next point-pair (step 440) . 
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When the user is finished, the system stores the active georeferencing 
functions with the raster-map (step is 445). At this time, the raster map is 
considered fully georeferenced. When accessed at any fiiture time, the system 
may simply retrieve the georeferencing functions, and apply them to fmd the 
latitude and longitude of any point on the raster map. 

The process of determining he georeferencing function set from a set 
of point-pairs is believed to be within the abiUty of one of ordinary skill in the 
art. The specific approach used by the system an method of the preferred 
embodiment is discussed below. 

At any given point in the georeferencing process a set of points? , has 
been specified. Let z e P denote one of the points of P . Associated value , 

, Latj and Lon. refer to thex and y coordinates of the point on the bitmap, 
and the latitude and longitude of the point on the earth's surface. We also 

define A c P to be the subset of "active" points, i.e. points which are 
currently being used in the calculation of the georeferencing functions. 
Assume that we have selected a parametric family of functions F. Functions in 
this family map from {x^y) to {Lon, Lat). Any function / eca/F has 
components fi{x,y) = Lon, and f2{x,y) = Lat . From within this family we 
seek the specific function / , which (in a sense to be defined below) comes 
closest to making the following system of equations true: 

f{^i>yi)'=^ {Lon,,Lat.) for i e A 

Once determined, / will b the georeferencing function which is used 
to compute corresponding latitude and longitude values, {Lon,Lat)for any 
point, on the bitmap. There are any number of possible ways to define 

the function that "comes closest to making (1) true." We shall follow a "least 
squares" approach also known in mathematics as an L^ approach. This 
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approach seeks to find the function, / , which minimizes the sum of the 
squared differences between the actual and the predicted values of latitude and 
longitude. In other words, firom among all the fimctions / e F, /is the one 
which minimizes: 

5 SSE = z[Mx,.yi)--Lon,f-i-(/^^^^^ (2) 

Among various alternative methods for choosing the function /are 
choosing it so t ht it minimizes the sum of absolute errors (rather than squared 
errors, or so that it minimizes the largest error. Other criteria are also possible. 

In what follows we shall describe methods for computing 
10 georeferencing functions based on both a "rotational linear" parametric family 

of functions and a "general linear" family of functions. A knowledgeable 
practitioner will readily perceive that this approach may be easily generaUzed 
to other parametric families of functions. Obvious candidates include 
nonlinear parametric famiUes of functions arising firom the projection of the 
1 5 earth' s surface onto a flat map. 

The General Linear Case: In the general linear case, we let F be the 
set of all possible linear transformations which map from (x^y) Xo{Lon,Lat) . 
Thus 



.^21 ^22 J W L4. 



(3) 



>V A 

20 for some choice of parameters , a^^^ a^x, ^22 » ^i- If the 

region covered by the map to be georeferenced is not too large, then this 
family of functions will contain a suitable function , / whose total error is 
quite small. In the case where the map to be georeferenced covers a larger area 
than this, then the curvature of the earth must be taken into account and F is 
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not a suitable family of fimctions. In such a case, nonlinear functions must be 
used as mentioned above. We shall not pursue that case further, since it is a 
straightforward extension of the procedures used in the linear case. 

To find / we seek the parameters which minimize 

5 SSE^Y^{a,,x, •^a.^y- -hb, -Lon.f +{a^,x. ^a^^y^ -hb^ -Lat.f The 

parameter values which minimize this expression are found by solving the 
following two independent systems of linear equations: 
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(5b) 



/eA 



These systems can be easily solve by well known methods, such as 
Gaussian EUmination, or LU factorization. The solutions yield the desired 

15 values of Oil , , ^21 > ^22* ^1 > ^2 • I* should be noted that equations (5a) 

and (5b) do not have a unique solution unless three or more non-colinear 
points are contained in A. Generally speaking, then, it requires 3 points to 
choose a georeferencing function from the family of general linear 
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transformations, When there are four points or more, it is possible to compute 
a standard deviation of errors using the formula: 



s = 



«-3 



(6) 



5 5 is an estimator for the amount of error to be expected between actual 

and predicted latitude and longitude values. 

Note that the inverse georeferencing function, /" that maps from 
(Lon, Lat) to (x,y) is readily obtained now, by inverting the function / . 
Having done this, it is possible to compute a similar standard deviation of 
10 error for which is an estimate of error as measured in the bitmap 

coordinates. 

The Rotational Linear Case: In the, so called, rotational linear case, 
we let F be the set of all possible linear transformations which map from 
{x,y)to {Lon, Lat) and which also allow for a known scale change caused by 
15 the difference in distance per degree of latitude and per degree of longitude at 

any point not on the equator. Thus 



-A 



r 

-A -A 



r 



(7) 



for some 



20 



choice of parameters >9i , A » A ' ^4 > where y = cos{Lat) is 
computed from the estimated latitude of the map to be georeferenced. 
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Clearly, this is a special case of the more general linear transform 
described above, where we take 

5 In theory, if the region covered by the map to be georeferenced is not 

too large, then this family of functions will contain a suitable function, 
/ whose total error is quite small. In the case where the map to be 
georeferenced covers a larger area than this, then the cxirvature of the 
earth must be taken into account and F is not a suitable family of 

10 functions. In such a case, nonlinear functions must be used as noted 

above. In practice, the general linear transform will yield results which 
are somewhat superior to the rotational linear transform due to the 
former's innate ability to use arbitrary scaling factors. These arbitrary 
scaling factors might compensate, for example, for scaling flaws in the 

15 digital scanner used to scan in the maps. The primary advantage of the 

rotational linear transform is that it can be computed with only two 
reference points (vmlike the general linear transform which requires 
three reference points). 

To find / we seek, as before, the parameters which minimize SSE as 
20 defined in equation (4). 

The parameter values which minimize this expression are foimd by 
solving the following system of linear equations: 
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These s>'stems can easily solved by well known methods, such as 
Gaussian Elimination, or LU factorization. The solutions yield the desire 
values of ^ h\y >9i which in turn yield the desired values for 

. , a^^, A, , and . 



ll should be noted that equation (8) does not have a unique solution 
10 unless two or m re points are contained in A. Generally speaking, then it 

requires two points to determine a georeferencing function from the family of 
rotational linear transformations. When there are three points or more, it is 
possible to compute a standard deviation of error, s using the formula: 



15 s is an estimator for the amount of error to be expected between actual 
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and predicted latitude and longitude values. 

Note that the inverse georeferencing function that maps from 

(Lon, Lat) to {xyy) is readily obtained now, by inverting the function /. 
Having done this, it is possible to compute a similar standard deviation of 
5 error for which is an estimate of error as measured in the bitmap 

coordinates. 

Automatic Error Detection and Handling When individual points 
are being assigned x, y, Lon, and Lat values, there is always a potential for 
error. To reduce the risk of incorrect georeferencing resulting from such 

10 errors, certain error handling procedures are built into the georeferencing 

process. The fundamental concem is that of detecting a "bad" point and then 
removing it from the set of active points, A Note that removing point from A 
will not delete the information associate, with that point, but it will cause the 
georeferencing parameters to be completely uninfluenced by that point. We o 

15 not wish to remove the point entirely, since it may be determined at a later 

stage of the georeferencing, that the point was not really bad at all, and should 
be-used in the georeferencing calculation. This will be clarified shortly. 

Detecting Bad Points The following steps outline the bad point 
20 detection process using the general linear transform approach to 

georeferencing. 

1 . Begin by placing all existing points into the active set, A 

2. If there are fewer than five active points then you are done. 
Otherwise, for each of the currently active points in tum, move it 

25 (call it point k for the sake of convenience) temporarily out of the 

active set, and then calculate the resulting inverse georeferencing 
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function (call it g^^*^) and its corresponding *SS£^jt : Also, calculate the 
difference between the predicted value and the actual value 

Sf^ =^g^^\Lonf.,Latf.)-{xf.^yf^)^. Make a note of the values, Jj^^^^*^^^/^^ • 

Return point k to the active set (and move on to the next value of k . 

5 3. From among the results found in step 2 above, find the 

3 / 

point, with the largest value of y^g^ )c^a?idS^ which also 

satisfies y^^^ yc^andSf^yc^^^here c^, andcj some constants which are 

set according to the general level of accuracy to be expected on the particular 
type of map which is being georeferenced, the current number of active points, 
10 and the dots per inch of the scanned image. If there is such a point then mark it 

as bad (by removing it firom the active set) and return to step 2 above. 
Otherwise you are done. 

There are several things to note about this procedure. One is that it 
allowing the value of c, and to change with the number of active points, 

15 makes it possible for the georeferencing system and method to utiUze points 

which it might originally determine bad or inconsistent after a large enough 
sample of points has been gathered to make it clear that a lesser level of 
accuracy is all that can be achieved on this map. Another observation is that 
using this procedure it is impossible to reduce the number of active points 

20 down to less than four (unless you started with less than 4 in which case this 

procedure does not apply at all). This scheme means that as each new point is 
added, all points determined so far are considered, even those which had 
previously been marked bad. Thus early "misjudgements" on the part of the 
system can be corrected later, in light of new point information, 

25 The same bad point detection process, can also be implemented using 

the rotational linear transform approach. In this case the method is capable of 
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reducing the niimber of active points down as low as three (rather than four for 
the general linear transform approach outlined above). This can be useful 
when dealing with small sets of active points. 

When the system is georeferencing raster map images that cover a 
large enough area so that the (nonlinear) ciuvature of the earth is a source of 
significant error, the system can calculate nonlinear georeferencing functions 
whose form corresponds to the map projection that was used to create Mapl. 

A specific example of the operation and application of the preferred 
georeferencing method may be shown with reference to the "Flood Zone 
Detertnination" business. The Federal Emergency Management Agency 
(FEMA) publishes a library of tens of thousands of paper maps showing 
various types of flood zones and their locations in the United States. A flood 
zone determination on a property is frequently done in the following way: 

1 . The address of the property is examined, and the 
location of the property is determined (perhaps through the 
use of a geocoding system, or b examining an available 
street map). 

2. A map analyst attempts to determine which of the 
many thousands of FEMA flood maps will contain this 
property. 

3. The map analyst goes to a map storage area and 
retrieves the desired map, often examining several maps 
before making a final selection. 
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4. Having retrieved the paper map, the map analyst 
next determines where, precisely, the property is located on 
the map. 

5. Finally, the map analyst examines flood zone 
notations on the map at the^rdpertys location in order to 
determine its flood-zone status. 

When performed using paper maps, the above process is difficult and 
quite time consuming . A database of scanned raster map images, alone, can 
be used to reduce the time and effort expended on step 3 above. However, 
georeferenced raster map images can also be used to eliminate virtually all 
work from steps 2 and 4 above. Since those are generally the most time 
consuming steps imder the current process, the value of georeferenced flood 
maps is considerable. 

Using georeferenced flood map raster images, steps 2 and 4 above, are 
replaced by: 

2. A computer system combines the pre-designated outlines of 
the raster map and the georeferencing information to obtain a polygon 
expressed in terms of latitude and longitude that outlines the 
region included in each flood map. Then the system determines which 
of the polygons contain the address in question, which is done using 
a "point-in-polygon" algorithm. At the conclusion of this process, the 
computer system has identified a map panel (or perhaps a small 
number of map panels) that contains the address. 
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4. Since the latitude and longitude of the property are known (by 
virtue of a geocoding phase), the computer system can use the 
georeferencing of the map panels to locate the property on each of the 
panels found above, thus largely eliminating any need for he map 
5 analyst to scan the flood map for the address location 

Therefore, it is easily seen that by using georeferenced raster maps 
instead of paper maps, the process of determining flood zones is substantially 
automated and much more efficient. 

It is important to note that while the present invention has been 
10 described in the context of a fully functional data processing system and/or 

network, those skilled in the art will appreciate that the mechanism of the 
present invention is capable of being distributed in the form of a computer 
usable medium of instructions in a variety of forms, and that the present 
invention applies equally regardless of the particular type of signal bearing 
15 medium used to actually carry out the distribution. Examples of computer 

usable mediums include: nonvolatile, hard-coded type mediums such as read 
only memories (ROMs) or erasable, electrically programmable read only 
memories (EEPROMs), recordable type mediums such as floppy disks, hard 
disk drives and CD-ROMs, an transmission type mediums such as digital and 
20 analog communication links. 

While the invention has been particularly shown and described with 
reference to a preferred embodiment, it will be understood by those skilled in 
the art that various changes in form and detail may be made therein without 
departing from the spirit and scope of the invention. 
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WHAT IS CLAIMED IS: 

1 . A method for georeferencing a raster map image, comprising 
the steps of: 

displaying a raster map and a georeferenced map; 
5 / identifying at least two geographically corresponding points on the 

raster map and on the georeferenced map; 

associating an image coordiDate of each point on the raster map with a 
geographic coordinate of the corresponding point on the goereferenced map; 
determining a functional relationship between the image coordinates 
10 and the geographic coordinates; and 

thereafter, for each additional corresponding points identified on the 
raster map and the georeferenced map, 

revising the functional relationship between the image coordinates and 
the geographic coordinates according to the additional corresponding points, 
15 and 

disregarding any points which are substantially inconsistent with the 
functional relationship. 

2. The method of claim 1, further comprising the step of: 
using the functional relationship to determine the geographic 

20 coordinates of features on the raster map. 

3. The method of claim 1, further comprising the step of: 
storing the functional relationship with the raster map. 

4. The method of claim 1, further comprising the step of: 
when the raster map is manipulated by a user, manipulating the 

25 georeferenced map accordingly. 

5. The method of claim 1, wherein the geographic coordinates are 
latitude and longitude. 
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6. The method of claim 1, wherein the raster map and the 
georeferenced map are displayed on the same computer display. 

7. The method of claim 1, wherein the corresponding points are 
marked by a user after visually detemiining geographically corresponding 
points. 

8. The method of claim 1 , wherein the functional relationship is 
represented by a set of general linear functions. 
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9. A computer system, having at least a processor connected to 
communicate with a readable and writeable memory, comprising: 

means for displaying a raster map and a georeferenced map; 

means for identifying at least two geographically corresponding .points 
5 on the raster map and on the georeferenced map; 

means for associating an image coordinate of the each point on the 
raster map with a geographic coordinate of the corresponding point on the 
georeferenced map; 

means for determining a functional relationship between the image 
10 coordinates and the geographic coordinates; and 

for each additional corresponding points identified on the raster map 
and the georeferenced map, 

means for revising the functional relationship between the image 
coordinates and the geographic coordinates according to the additional 
15 corresponding points, and 

means for disregarding any points which are substantially inconsistent 
with the functional relationship. 

10. The system of claim 9, further comprising: 

means for using the functional relationship to determine the geographic 
20 coordinates of features on the raster map. 

1 1 . The system of claim 9, further comprising: 

means for storing the functional relationship with the raster map. 

12. The system of claim 9, further comprising: 

means for, when the raster map is manipulated by a user, manipulating 
25 the georeferenced map accordingly. 

13. The system of claim 9, wherein the geographic coordinates are 
latitude and longitude. 
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14. The system of claim 9, wherein the raster map and the 
georeferenced map are displayed on the same computer display. 

15. The system of claim 9, wherein the corresponding points are 
marked by a user after visually determining geographically corresponding 

5 points. 

16. The system of claim 9; wherein the functional relationship is 
represented by a set of general linear functions. 
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17. A method of georeferencing a raster map, comprising: 
receiving a display of a first map and a second niap, the first map being 

a digital raster map, the second map being a previously georeferenced map; 

receiving a manipulation of the first map such that the first map and 
the second map are approximately aligned; 

displaying the first map and the second map; 

receiving a selection of a point pair point on the first map; and 

receiving a selection of a corresponding point pair point on the second 

map. 

1 8. The method of claim 1 wherein the second map is a vector 

map. 

1 9. The method of claim 1 wherein the second map is a digital 
raster map. 

20. The method of claim 1 wherein the point pair point on the first 
map has a known longitude and latitude. 

21. The method of claim 1 further comprising creating a 
georeferencing function. 

22. The method of claim 5 wherein the georeferencing fimction is a 
linear transformation. 

23. The method of claim 1 further comprising selectively 
synchronizing, responsive to a user command, the first map and the second 
map. 

24. The method of claim 1 fiirther comprising receiving a mark on 
the first map at a location, and reproducing the mark on the second map at a 
corresponding location. 



wo 01/73690 



PCT/USOl/10154 



25. The method of claim 5 wherein georeferencing uses at least 
three point pairs to complete the georeferencing function for the first map 
based on a linear transformation. 

26. The method of claim 5 further comprising using at least four 

5 point pairs to complete the georeferencing ftmction for the first map, based on 

a linear transformation, and fiirther comprising executing a validation check. 

27. The method of claim 10 further comprising rejecting a point 
pair u hen the point pair deviates a pre-determined amount fi"om a pre- 
dclcmuncJ standard error. 

10 28. The method of claim 1 1 wherein the predetermined standard 

error uses a "least square" parameter fitting operation. 

29. The method of claim 1 fiirther comprising : 
receiving a selection of a point pair on the first map, and 
receiving a selection of a point pair on the second map. 
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30. A method in a computer system for georeferencing a raster 
map, comprising: 

providing for display a first map and a second map, the first map being 
a digital raster map, and the second map being a previously georeferenced 
map; 

the fist map being substantially similar to the second map when 
displayed; 

receiving an entry identifying a first point pair point on the first map; 

receiving an entry identifying a second point pair point on the second 
map, the second point pair point having approximately the same location on 
the second map as the first poiijt pair point has on the first map; 

assigning a point pair point on the first map a longitude coordinate and 
a latitude coordinate, the longitude coordinate and the latitude coordinate of 
the first point pair point being identical to a longitude point and a latitude 
point associated with a point pair point on the second map. 

3 1 . The method of claim 1 4 further comprising creating the 
georeferencing function, 

32. The method of claim 14 further comprising receiving a mark on 
the first map at a location, and reproducing the mark on the second map at a 
corresponding location. 

33. The method of claim 14 further comprising using at least four 
point pairs to compute a georeferencing function for the first map based on a 
linear transformation, and further comprising executing a validation check. 

34. The method of claim 17 further comprising rejecting a point 
pair when the point pair deviates a predetermined amount firom a 
predetermined standard error. 
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35. A computer readable medium whose contents enable the 
georeferencing of a raster map, by: 

providing for display a first map and a second map, the first map being 
a digital raster map, and the second map being a previously georeferenced '■ 
5 map; 

the fist riiap being substantially similar to the second map when 
displayed; 

receiving an entry identifying a first point pair point on the first map; 
receiving an entry identifying a second point pair point on the second 
10 map, the second point piair point having approximately the same location on 

the second map as the first point pair point has on the first map; 

assigning a point pair point on the first map a longitude coordinate and 
a latitude coordinate, the longitude coordinate and the latitude coordinate of 
the first point pair point being identical to a longitude point and a latitude 
15 point associated with a point pair point on the second map. 

36. The method of claim 19 fiirther comprising: 

using at least four point pairs to compute a georeferencing fimction for 
the first map based on a hnear transformation; 

fiirther comprising executing a validation check; and 
20 rej ecting a point pair when the point pair deviates a predetermined 

amount fi-om a predetermined standard error. 
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37. A method of georeferencing a raster map, comprising; 
receiving a display of a first map and a second map, the first map being 
a digital raster map, the second map being a previously georeferenced map; 
receiving, a manipulation of the first map such that the first map and 
5 the second map are approximately aligned; 

displa>'ing the first map and the second map; 

receiving a selection of a point pair point on the first map; and 

rcx^eivmg a selection of a corresponding point pair point on the second 

map. 

10 38. The method of claim 1 fiirther comprising verifying that the 

point pair point on the first map is correctly associated with the point pair 
point on the second map. 

39. The method of claim 1 wherein aligning places the first map 
within the second map. 

15 40. The method of claim 1 fiirther comprising the act of providing 

a longitude and latitude to the point pair point on the second map. 

4 1 . The method of claim 1 wherein the point pair point on the 
second map has a known longitude and latitude. 

42. The method of claim 1 fiirther comprising generating a 
20 georeferencing function. 

43. The method of claim 1 fiulher comprising providing an 
identification of a plurality of point pair points. 

44. The method of claim 1 fiirther comprising the act of marking a 
point on the first map, the point on the first map being automatically 

25 reproduced on the second map. 
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45. The method of claim 8 further comprising correcting the 
reproduced mark. 

46. The method of claim 1 further comprising selecting a 
predefined georeferencing function to associate a point on the first map with a 

5 point on the second map. 
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47. A computer readable medium whose contents transform a 
computer system into a raster map georeferencing system, by enabling: 

receiving a display of a first map and a second map, the first map being 
a digital raster map, the second map being a previously georeferenced map; 

manipulating the first map such that the first map and the second map 
are approximately aUgned; 

selecting a point pair point on the first map; and 

selecting a corresponding point pair point on the second map. 

48. The computer readable medium of claim 1 1 wherein the 
contents of the computer readable medium are also capable of verifying that 
the point pair point on the first map is correctly associated with the point pair 
point on the second map. 

49. The computer readable medimn of claim 1 1 , wherein the 
contents of the computer readable medimn are also capable of incorporating a 
plurality of point pair points. 

50. The computer readable medium of claim 1 1 wherein the 
contents of the computer readable medium are also capable of allowing a user 
to mark a point on the first map, the point on the first map being automatically 
reproduced on the second map. 

51. The computer readable medium of claim 1 1 wherein the 
contents of the computer readable medium are also capable of providing a 
longitude and latitude to the point pair point on the second map. 
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52. A computer memory containing a data structiire capable of 
enabling the georeferencing of a raster map, by: 

receiving a display of a first map and a second map, the first map being 
a digital raster map, the second map being a previously georeferenced map; 
5 receiving a manipulation of the first map such that the first map and 

the second map are approximately aligned; 

displaying the first map and the second map; 

receiving a selection of a point pair point on the first map; and 

receiving a selection of a corresponding point pair point on the second 

10 map. 

53. The computer memory of claim 16 wherein the computer 
memory maintains code that verifies that the point pair point on the first map 
is correctly associated with the point pair point on the second map. 

54. The computer memory of claim 16, wherein the computer 
15 maintains code that supports the providing of a parality of point pair points. 

55. The computer memory of claim 16 wherein the computer 
memory maintains code for marking a point on the fiirst map, the point on the 
first map being automatically reproduced on the second map. 

56. The computer memory of claim 16 wherein the computer 

20 memory maintains code that provides longitude and latitude to the point pair 

point on the second map 
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57. A system that enables the georeferencing of a digital raster 
map, comprising: 

a processing platform for executing code capable of georeferencing a 
digital raster map; and 
5 a storage platform for storing at least a digital raster map, the storage 

platform coupled to the processing platform. 

58. The system of claim 1 further comprising a user interaction 
device coupled to the processing platform. 

59^ The system of claim 1 wherein the processing platform is 
10 enabled to georeference a digital raster map based on a second map that is 

already georeferenced. 

60. The system of claim 1 wherein the processing platform is a 
microprocessor. 

61 . The system of claim 1 wherein the map processing platform is 
15 a application service provider. 

62. The system of claim 1 wherein the map processing platform is 
located remotely from a user of the processing platform. 

63. The system of claim 1 wherein the storage platform comprises 
cashed memory. 

20 64. The system of claim 1 wherein the storage platform comprises 

system memory. 

65. The system of claim 1 wherein the storage platform comprises 
non-cashed volatile storage. 
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66. The system of claim 2 wherein the user interaction device 
comprises at least a display device. 

67. The system of claim 1 where in the processing platform is 
coupled to the storage platform via a network. 

5 68. The system of claim 1 1 wherein the network is the intemet. 

69. The system of claim 1 wherein the storage platform maintains a 
parality of digital raster maps. 

70. The system of claim 13 wherein the storage platform maintains 
a database of georeferencing functions that associate a digital raster map 

10 location with a georeferenced location in a second map. 

71 . The system of claim 1 wherein the storage platform maintains 
code that enables the georeferencing of a raster map, by: 

receiving a display of a first map and a second map, the first map being 
a digital raster map, the second map being a previously georeferenced map; 
15 receiving a manipulation of the first map such that the first map and 

the second map are approximately aUgned; 

displaying the first map and the second map; 

receiving a selection of a point pair point on the first map; and 

receiving a selection of a corresponding point pair point on the second 

20 map. 
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72. A data signal comprising a data structure capable of 
georeferencing a raster map, by: 

providing for display a first map and a second map, the first map being 
a digital raster map, and the second map being a previously georeferenced 
5 map; and 

assigning a point pair point on the first map a longitude coordinate and 
a latitude coordinate of a point pair point on the second map. 

73. The method of claim 16 further comprising creating a 
georeferencing function. 

1 0 74. The method of claim 1 6 further comprising receiving a mark on 

the first map at a location and reproducing the mark on the second map at a 
corresponding location. 

75. The method of claim 16 further comprising using at least four 
point pairs to complete the georeferencing function for the map based on a 

15 linear transformation, and further comprising executing a validation check. 

76. The method of claim 1 9 further comprising rejecting a point 
pair when the point pair deviates a pre-determined amount firom a pre- 
determined standard error. 
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